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Abstract 

Using a generating function approach, the correlation coefficients of 
four different graph-theoretical indices, namely the number of independent 
vertex subsets, the number of matchings, the number of subtrees and the 
Wiener index, are asymptotically determined for random rooted ordered 
trees. 



1 Introduction 

In the author investigates correlation measures for graph-theoretical indices 
which are of interest in theoretical chemistry. In particular, the correlation 
coefficients for these indices are asymptotically determined. Since the necessary 
calculations are rather lengthy and tedious, only a few of them are explicitly 
provided there. This additional note fills the gap by explaining the involved 
details. See m] for further references and applications. 

The underlying stochastic model is the following: of all rooted ordered trees 
on n vertices, a tree T n is selected uniformly at random. The parameter we 
are interested in is the correlation coefficient of two indices X n — X(T n ) and 
Y n = Y(T n ), defined by 

v/Var(X„)Var(r„) 
Here, X and Y are two of the following four indices: 

(1) The Merrifield- Simmons- or a- index is defined to be the number of inde- 
pendent vertex subsets of a graph, i.e. the number of vertex subsets in 
which no two vertices are adjacent, including the empty set. Merrifield 
and Simmons investigated the er-index in their work [5] and pointed out 
its correlation to boiling points of molecules. 
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(2) The Hosoya- or Z- index is denned as the number of independent edge 
subsets (also referred to as "matchings" ) , i.e. the number of edge subsets 
in which no two edges are adjacent, including the empty set again. 

(3) The number of subtrees is called p-index in [8] and was discussed lately in 
a paper of Szekely and Wang |1U) . 

(4) The Wiener index is probably the most popular topological index (s. 
\V2\). It is defined as the sum of all the distances between pairs of vertices, 
i.e. 

W(G)= J2 d oM- (2) 

v,w<£V{G) 

It will be shown how one can derive functional equations for generating 
functions yielding to asymptotic formulas for the quantities of interest by means 
of the Flajolet-Odlyzko singularity analysis 3 . It turns out that the approach 
is a slightly different one for the correlation with the Wiener index in view of 
its different growth structure. Therefore, the correlation coefficients of the a-, 
Z- and p-index are determined in the following section, whereas the correlation 
to the Wiener index is investigated in Section [3] 

2 a-, Z- and p- index 

In this section, we want to determine the asymptotic behavior of the generating 
function 

T 

where X, Y stand for a-, Z- or p-index (possibly, X and Y are the same). They 
count, respectively, the number of independent vertex subsets, edge subsets and 
subtrees of a tree T. This is done by distinguishing between two cases for each 
of the indices: 

• the root vertex belongs to the independent vertex subset / edge subset / subtree, 

• the root does not belong to it. 

We denote the corresponding quantities by o\ (T) , a% (T) resp. Z\(T), Z 2 {T) and 
pi(T), p2(T). If T%, . . . ,Tk are the branches of the rooted tree T, it is easy to 
see that 

k 

(T 1 (T) = l[a 2 (T i ), 

k 

o-2{T) ^Xlio-^) + o^T,)), 

i=l 



2 



k 

Z 2 (T) = l[(Z 1 (T i ) + Z 2 (Ti)), 

k 

i=l 
k 



i-s 2 ( z y 

z 



The corresponding generating functions are called S\,S 2 resp. Zi,Z 2 and 
R\ , R 2 ■ Functional equations for these functions which follow from the recursive 
relations given above and lead to asymptotic formulas for the average indices 
have already been presented by Klazar [7] and others. For the sake of complete- 
ness, these are given here as well: 

Si(z 
S 2 (z 

Zi(z 
Z 2 (z 

Ri(z 

R 2 {z 



1 - 


Si(z) - S 2 {z) ' 




zZ 2 (z) 


(1- 


-Z l {z)-Z 2 {z)f 




z 


1 - 


z 1 { z )-z 2 { z y 




z 



l-R 1 {z)-T{ Z y 



:(Ri(z) + R 2 (z)). 



Here, T(z) is the generating function for the number of rooted ordered trees. It 
is well-known that T(z) satisfies the functional equation 

T(z) 



i-T( z y 

which leads to an explicit formula for the number of rooted ordered trees: 



t n = 



1 (2n - 2 



n V n — 1 



1 



40? 



-3/ 2/ in 



As determined by Klazar , the functional equations presented above yield the 
following asymptotic formulas for the expected values of our indices: 



E{a n ) ~ \/3 



(1.02640) • (1.6875)™. 
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E(Z n ) 



'65 



78 



13 / 35 + 13 V / 13 N 



54 



(0.88719) • (1.51615)™. 



E(Pn) 



16 



3V15 



— w (1.37706) • (1.5625) 1 



If we combine cr^a and ^1,2, we obtain generating functions SZn, SZ12, SZ21, SZ22 
for the correlation between a- and Z-index: 



SZij(z) 



T 



<j t {T)Z 3 {T)z 



\T\ 



The total is denoted by SZ. In an analogous manner, we define SR^ ,ZRjj ,SSij ,ZZj 
and RRjj . 

Next, we observe that 1 < a(T), Z(T), p{T) < 2^ for trivial reasons. This 
helps us to restrict the range of the radius of convergence. In all our cases, it 
must lie within the interval [yg , |] . Even more can be told about it: the radius of 
convergence of XY can be at most the minimum of the radii of X and Y (which 
are ^, 13v ^f~ 35 and ^ for a-, Z- and p-index respectively). Furthermore, 
since <r(T) > F\t\+2, where F n denotes the n-th Fibonacci number (a result 
due to Prodinger and Tichy [Hj), we even know that, for example, the radius of 



VE-l _ 2(x/5-l) 



For the generating functions SS, 



convergence of SR is < 9 . , — .,. 
ZZ and RR, we may apply the Cauchy-Schwarz-inequality to obtain an upper 
bound for the convergence radius easily: if X has a convergence radius r, then 
XX has convergence radius < 4r 2 . Summing up, we have the following bounds 
for the convergence radii: 



SZ: interval 



(13-y/l3-35)(v / 5-l) 



144 



J_ 2(V5-1) 
16' 25 



• SR: interval 

• ZR: interval [^,^], 

• SS: interval [^3, 7^], 

• ZZ: interval 



j_ 1711-455VT3 

16 ' 648 



• RR: interval [j^pggg]- 

In all the cases, we will see that these estimates are sufficient to determine the 
correct dominating singularity. 
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2.1 a- and Z-index 

The recursive relations for a- and Z-index lead to the following system of func- 
tional equations: 

zSZ 22 (z) 



SZ u (z) 
SZi 2 (z) 

SZ 21 (z) 
SZ 22 (z) 



(l-SZ 21 (z)-SZ 22 (z)) 2 ' 

z 

l-SZ 21 (z)-SZ 22 (z)' 

z(SZ 12 (z) + SZ 22 (z)) 

(1 - SZ n (z) - SZ 12 (z) - SZ 21 (z) - SZ 22 (z)) 2 ' 

z 



1 - SZn(z) - SZi 2 (z) - SZ 21 (z) - SZ 22 (z) 
For instance, the functional equation for SZn can be deduced as follows: 

SZ 1 i(*)=$>i(r)Zi(:Z> |T| 



T 

ft 



2 |T!|+...+|T fe |+l 



^ k SZ 22 (z)(SZ 21 (z) + SZ 22 ( Z )) fc - 1 



z 

k>0 

zSZ 22 (z) 



(l-5Z 21 (^)-SZ 22 (z)) 2 - 

From these, a single equation for SZ 22 can be worked out by means of Grobncr 
bases 0], and this can be used to determine the dominating singularity of SZ. All 
computations are given in the accompanying Mathematica® files, which can be 
found on http://finanz.math.tugraz.at/~wagner/Correlation. From the 
equation 

s 10 + 2zs 8 - 3zs 7 + z 2 s 6 - 4z 2 s 5 + 3z 2 s 4 - z 3 s 3 + 2z 3 s 2 - z 3 s + z 4 = 
that is satisfied by s = SZ 22 (z), we find that 



SZ 22 (z) - 0.171502 - 0.138532. /l - — 

V z o 

around the singularity zq = 0.0982673. From the relation SZ(z) = 1 — g^Tz) 
we obtain 

SZ0) - 0.427020 - 0.462827. Il- — , 

V z o 

which gives us the asymptotic formula 

E(a n Z n ) ~ (0.92565) • (2.54408)". 

by a simple application of the Flajolet-Odlyzko singularity analysis as it is 
also explained in |11|. 
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2.2 a- and p-index 

The recursive relations for a- and p-index lead to the following system of func- 
tional equations: 

SRn(z) 

SRi 2 (z) 
SRai(«) 
SR 22 (z) 



l-SR 21 (z)-S 2 (z)' 

^(SR 21 (z) + SR 22 (z)) ^(z) 2 (SR 21 (z) + SR 22 (z)) 
(I-S2W) 2 

z 

1 - SRn(z) - SR 21 (z) - 5i(z) - S 2 (z)' 
z(SR n (z) + SR 12 (z) + SR 21 (z) + SR 22 (z)) 

(1 - S 1 (z) - S 2 (z)) 2 
S 2 (z) 2 (SRn(z) + SR 12 (z) + SR 21 (z) + SR 22 (z)) 



It would be possible to carry out the same procedure as in the previous case; 
however, it saves computing time to consider SRn and SR 2 i first. Then, SRi 2 
and SR 22 are given by simple linear equations which result in the formula 

m , z(S 1 (z) 2 SR 21 (z) + z SR n (z) + z SR 21 (z)) 
(Z) z 2 - zS 2 {z) 2 - S 1 {z) 2 S 2 {z) 2 

We know that S\ and S 2 are holomorphic in the region of interest, so the 
dominating singularity of SR is either a singularity of SRn and SR 2 i or a zero 
of the denominator. However, using the functional equations for Si and £2, we 
find that the denominator only vanishes at z = and at z = Since ^ does 
not lie in our estimated interval, we have to determine the singularities of SRn 
and SR21. By the Grobner basis approach, we find that s = SRn(z) satisfies 
the polynomial equation 

zs 9 - 6zV - 4z 3 s 6 + 7z 3 s 5 - 2z 4 s 4 - (z 5 + 3z 4 )s 3 + z 5 s 2 + z 5 s - z 6 = 0. 

is clearly not a singularity, and the other common zeroes of the polynomial 
and its derivative satisfy the polynomial equation 

(27z - 4)(8z 2 + 81) 2 (125z 3 - 412z 2 - 40936z + 3844) = 0. 

The only value that lies within our bounds and satisfies the equation is z ~ 
0.0938166. Since SR 2 i(z) = 1 — S 2 (z) — SR ^ ^ and SRn only vanishes at z = 0, 
the smallest singularity of SR 2 i is the same. Expanding around z gives us 



SR(z) - 0.560623 - 0.683264. /l - — 

V z o 

and thus 

E(a nPn ) ~ (1.36653) • (2.66477)™. 
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2.3 Z- and p-index 

The recursive relations for Z- and p-index lead to the following system of func- 
tional equations: 



ZRn(z) 
ZR 12 (z) 

ZRai(z) 
ZR 22 (z) 



z(ZR 21 (z) + Z 2 (z)) 



(1 - ZRn(z) - ZR 21 (z) - Z x {z) - Z 2 {z)) 2 ' 
z(ZR 2i (z) + ZR 22 (z)) 2zZ 2 (z)ZR(z) 



(1-Z!(Z)-Z 2 (Z)) 2 (1-Zl(*)-Z 2 (*)) 3 

Z 2 (z) 2 (ZR 21 (z) + ZR 22 (z)) 2Z!(z)Z 2 (z)ZR(z) 



1 - ZRn(z) - ZR 2 i(z) - Z x {z) - Z 2 (z) ' 
zZR(z) Z 2 (z) 2 ZR(z) 



(1-Z 1 (z)-Z 2 (z)) 2 z 
Using the same approach as in the previous example, we arrive at 
z{Z 2 {zf ZR 2 i(z) + z ZRn(z) + zZR 21 (z)) 



ZR(z) 



z 2 - 2zZ 1 (z)Z 2 (z) - zZ 2 (z) 2 - Z 2 (z) 4 



Again, Z\ and Z 2 are holomorphic in the region of interest, and the denominator 
only vanishes at z — and at z = 13v/ I|~ 35 , which does not lie in our estimated 
interval, so we have to determine the singularities of ZR n and ZR 21 . We observe 
first that ZR n (z) = ZR2i(*) 2 (ZR2i(*)+z 2 (*)) an d Z x {z) = ^1. Inserting yields 

z 2 - z ZR 21 (z) + zZ 2 (z) ZR 2 i(z) + Z 2 {zf ZR 21 (z) 

+ zZR 21 (z) 2 + Z 2 (z) ZR 21 (z) 3 + ZR 21 (z) 4 = 0. 

Elimination of Z 2 via the relation z 2 — zZ 2 (z) + zZ 2 (z) 2 + Z 2 (z) A = gives 
a polynomial equation of degree 16. For computational purposes, however, it 
is much faster to find a common zero of the equation above together with its 
derivative and the condition for Z 2 . Then we find that a singularity of ZR 2 i 
(and thus also of ZRn) and ZR must satisfy the polynomial equation 

(4096z 2 - 448z + l)(2560000z 2 + 2894400z + 531441) = 0. 

The only value that lies within our bounds and satisfies the equation is zo ~ 
0.107095. Expanding around z = gives us 



zr W~t4( 211 + 93 ^)-^ 



'5(128985 + 57683^5) 



1-^ 



928 v 7 232 V 58 V z o 

and thus 
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2.4 Variance of the cr-index 

For the variances, the calculations arc even a little simpler, since we have one 
variable less to cope with. In this case, we obtain the functional equations 

SSii{z) = i-ssm> 

SS 12 (z) 



1 - SS 12 (z) - SS 22 (z) ' 

SS 22 (z) 



1 - SSu(z) - 2SS 12 (z) - SS 22 {z) ' 
which result in a single equation for s — SS(z) = SSn(z) + 2SS\ 2 {z) + SS 22 (z): 

s 6 - 6s 5 + (4z + 14)s 4 + (8z 2 - 20z - 16).s 3 + (4z 3 - 30z 2 + 36z + 9)s 2 

- (12z 3 - 36z 2 + 28z + 2)s - (z 4 - 8z 3 + 14z 2 - 8z) = 0. 

We use Grobner bases once again and find the only possible singularity which 
lies within our bounds: its value is z ~ 0.0873832. We calculate the expansion 
of SS(z) around z : 



SS{z) ~ 0.614803- 0.519010. /l - — , 

V z o 

and finally obtain the asymptotics of the variance: 

Var(<T„) - (1.03802) • (2.86096)". 

2.5 Variance of the Z-index 

We proceed in the same way in the case of the Z-index: 

zZZ 22 (z) 2z(ZZ 12 (z) + ZZ 22 (z)) 2 



ZZ n (z) = 
ZZ 12 (z) = 



(l-ZZ(s)) 2 (l-ZZ(s)) 3 
z(ZZ 12 (z) + ZZ 22 (z)) 



(l-ZZ(z))2 ' 

zz - (z) = T=mry 

which results in a single equation for s = ZZ(z) = ZZ 1:L (z) + 2ZZ 12 (z) + ZZ 22 (z): 

s 8 - 7s 7 -(z- 21)s 6 + (4z - 35)s 5 + (2z 2 - 5z + 35)s 4 - (7z 2 + 21)s 3 

- (z 3 - 9z 2 -5z- 7)s 2 + (2z 3 - 5z 2 - 4z - l)s + (z 4 - z 3 + z 2 + z) = 0. 

Here, the value of the singularity is z ~ 0.107969. We calculate the expansion 
of ZZ(z) around z : 



ZZ(z) - 0.296221 - 0.386136 Jl - — , 

V z o 

and finally obtain the asymptotics of the variance: 

Var(Z„) - (0.77227) • (2.31549)™. 
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2.6 Variance of the p-index 

In this case, we obtain the following system of equations: 

z 



RRn(z) 
RR 12 (z) 



(1 - RRn (z) - 2i?i (z) - TO) ' 
z(RRn(z) + RRigjf) + Ri(z) + R 2 (z)) 

(l-R^-TizW 
i?i(z) 2 (RR n (z) + RRi 2 (z) + R^z) + R 2 {z)) 



z 

^_ zRR(z) 2zQR 1 (z) + fl 2 (z)) 2 

RR22(z) - d-T{z)f + — (i-T(z)r — 

_ T(z) 2 RR(z) | 2i? 2 (z) 2 (l - T(z)) 

z ' z 

Note that the system can even be solved explicitly by successive solution of 
quadratic equations. We apply the method that was also used for the covariance 
of a- and p- resp. Z- and p-indcx: RR i2 and RR 22 are expressed in terms of 
the other functions. We obtain 

N 

RR(z) = RRn(z) + 2RR 12 (z) + RR 22 (z) = — 2 z)(^(z) 2 - z) ' 

where the numerator N is given by 

N = zRR n (z)(z + i?i(z) 2 ) - 2i? 1 (z) 2 i? 2 (z) 2 (l - T(z)) 

+ 2zi? 1 (z) 2 (i? 1 (z) + R 2 (z)) + 2zR 2 {zf{\ - T(z)). 

Next, we prove that the denominator doesn't vanish within the bounds of in- 
terest: an easy calculation shows that it can only vanish at z = 0, z = \ or 
z — — 

Z ~ 25- 

So it suffices to determine the singularities of RRn(z). The application of 
Grobner bases shows that the smallest singularity of RRn(z) is the surprisingly 
nice value zq = Expansion of RRn around zq yields 

RRu(,)~^- * /TZ 
V ; 9 3V7V V 

so after calculating the values of Ri, R 2 and T at z , we see that 



2432V2 - 1632 32\/l4 

giving us the asymptotics 



64VT4 {81 

Var(p„) 



147 V 32 
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3 The Wiener index 

Now, we want to determine the asymptotic behavior of the generating function 

^X(T)^(T)zl T l, 

T 

where X stands for a-, Z- or /9-index,VF(T) is the Wiener index and T runs 
over all rooted ordered trees. Again, we have to distinguish sets containing and 
not containing the root. For the Wiener index, on the other hand, we have the 
following recursive relations: 

k 

D(T) = J2D(Ti) + \T\-l 

i=i 

and 

k 

W{T) = D(T) + W(Ti) + (p(Ti) + |T,|) \Tj\. 

Here, D is the total height or internal path length, the sum of the distances to 
the root. From these relations, one finds the following functional equations for 
the respective generating functions: 



zD{z) 



D{z) = - ^ > +zT'(z)-T(z), 



w{z) = D(z) + (1 _ r(z))2 + ( i-r(z))3 ' 

from which D(z) and W(z) can be determined without difficulty. The asymp- 
totic behavior of the average Wiener index follows at once: 

E(W n ) ~ 

The variance of the Wiener index is given in a paper of Janson [S]: 

16- 5tt 



Var(VF„) 



80 



-n 



5 



Now, we define generating functions DSi, DS2, WSi, WS2 for the product of 
D(T) resp. W(T) with the number of independent vertex subsets containing 
resp. not containing the root, e.g. 



DS 1 (z) = J2D(T)a 1 (T) 



In an analogous manner, we define the functions DZ^ , WZ^ , DR^ , WRi . Here, 
we obtain linear functional equations for the generating functions, which can be 
solved explicitly. 
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3.1 a- and Wiener index 

The recursive relations give us the following system of functional equations, 
which can be simplified by means of the functional equations for S± and S 2 , 
especially the facts that 1 _| 2(z) = Si(z), 1 _ Sl l S2(z) = S 2 (z) and (1 _g^ (z))2 = 
52(2) (the latter follows after some simple algebraic manipulations): 



DS 2 (z)S 2 (z) + zS[(z) - 5i(z) 

z(DSi(z)+DS 2 (z)) 
(1-S 1 (z)-S 2 (z)) 

S 2 (s) 2 (DSi(s) + DS 2 (*)) 



DS 2 (.) = + *S' 2 {z) - S 2 (z) 



+ zS' 2 (z) - S 2 (z), 

ww zWS 2 (z) 2z*S' 2 {z){m 2 {z) + zS' 2 {z)) 
W Sl (z) = D Sl (z) + (1 _ Sa(z))2 + (1-S 2 (z))3 

= DSi(«) + S 2 (*) WS 2 (z) + 2^(z)5 2 (z)5 2 (z)(DS 2 (z) + z5 2 (z)), 

z(WSi(z)+WS 2 (z)) 



WS 2 (z) =DS 2 (z 
+ 



(1-Si(z)-S 2 (*)) 2 
2z(^(z) + zS 2 (z))(DSi(z) + DS 2 (z) + z^(z) + zS' 2 {z)) 



= DS 2 («) 



(1-^(^-52(^)3 

5 2 (z) 2 (WS!(z)+WS 2 (z)) 



z 

| 2S 2 (z) 3 (^(z) + S 2 (z))(DS!(z) + DS 2 (z) + zS^(z) + zS' 2 {z)) 

z 

All these equations are obtained as in the following example: 



|T| 



DSi(z) = E £) ( I >i( 7 > 



= EE-E E II ^) - |Til+ - +|Tfcl+1 
fe>o Ti T fc \i=i i=i y 

+ ^(|T|-l)a 1 (T)^l 

T 

= EE-E fe^m^mjii^w)) - |Til+ - +|Tfcl+1 

fe>o Ti T fc \i=i y 
+ *S , 1 (*)-S , i(z) 
= «^ fcDS 2 (z)5 2 (z) fc - 1 + zS[(z) - S^z) 

k>0 

= {i-s 2 My +zs[{z) - si{z) - 
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We solve the system of linear equations for WSi and WS2 and obtain an expres- 
sion for WS(z) = WSi(z) + WS 2 (z) in terms of Si and S 2 . Then, we replace 
^by^and S[(z) by 

S[{Z) = ll-S 2 (z) = l-l 2 (z) + (l^tiz)? = l-S 2 {z) + S ^ S '^- 

Finally, implicit differentiation of the equation S 2 (z) 3 — 2S 2 (z) 2 + S 2 (z) — z = 
yields 

S ' 2(Z) = 35 2 (z) 2 -45 2 (z) + l' 

so WS can be written in terms of £2 and z only. We obtain an expression of 
the form 

N 

WS(z) = 



(1 - 3S 2 (z)) 2 (l - S 2 {z)f(S 2 {z) 2 + S 2 (zf - zf ' 



where N is a polynomial in S 2 and z. The denominator only vanishes at z = 
(which is clearly no singularity) and at z = which is the dominating 
singularity of S 2 . Therefore, we only have to consider the expansion of WS 
around which is given by 

WS{z) 



2 ■ 



81(1-^) 
This yields 

20V5F 5/2 (27 



3.2 Z- and Wiener index 

All steps are analogous to the previous section. We start from the following func- 
tional equations, which arc simplified by the relations Z\(z) = (i^z 1 (z)-z 2 (z))' i 
and Z 2 (z) = ^z^-z^y 



D7 M 2zZ 2 (z)(T>Z 1 (z)+'DZ 2 (z)) zDZ 2 (z) 

DZ1(Z) - (1 - Z l{ z) - Z 2 {z)f + (1-Z l{ z)-Z 2 (z)) 2 + ZZ1(Z) ~ Z1(Z) 



2Z 1 (z)Z 2 (z)(DZ 1 (z) + DZ 2 (z)) Z 2 (z) 2 DZ 2 (z) 

H h zZ 1 (z) - Zi(z), 



z z 



z(DZ 1 (z) + DZ 2 (z)) 

DZ2(z) - F^FW + 2( } " 2( } 

^ 2 (,) 2 (DZ l( ,) + DZ 2 (,)) + ^ (z) _ Z2(z); 
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WZ l( ,) = DZ l( ,) + 2 ^ W(WZlW+WZ2(2)) ' ZWZ2(Z) 



+ 



+ 



(1 - Z 1 (z) - Z 2 (z)f (I - Z^z) - Z 2 {z)Y 

Or / 

■ ((DZ 2 (z) + zZ' 2 (z))(zZ[(z) + zZ' 2 (z)) 



(1-Z 1 (z)-Z 2 (z))3 

+ zZ 2 (z)(DZi(z) + DZ 2 (z) + zZ[(z) + zZ' 2 {z)y 

6zZ 2 (z)(zZ[(z) + zZ£(z))(DZi(z) + DZ 2 (z) + zZj(z) + zZ 2 (z)) 
(1-Z 1 (z)-Z 2 (z))* 



= DZi(z) + 2Z 1 (z)Z 2 (z)(WZ 1 (z)+WZ 2 (z)) + Z 2 (z) 2 WZ 2 (z) 

z z 
+ 2Zx(z)((DZ 2 (z) + zZ 2 (z))(Z[(z) + Z' 2 {z)) 

+ Z 2 (z)(DZ!(z) + DZ 2 (z) + zZ[(z) + zZ' 2 {z))) 
+ QZl(z) ^ {z) \ z[{z) + Z' 2 {z))<J>Mz) + DZ 2 (z) + zZ[(z) + zZ' 2 {z)), 
WZ 2 (z)^DZ 2 (z)+ Z ( WZl W +WZ2 W) 



+ 



{l- Zl {z)-Z 2 {z)f 
2z 2 {Z' l {z) + Z£(z))(DZi(*) + DZ 2 (z) + zZ' x (z) + zZ' 2 (z)) 



DZ 2 (z) + 



{\-Zy{z)-Z 2 {z)f 

Z 2 (z) 2 (WZ 1 (z)+WZ 2 (z)) 



2Z 2 {zf{Z[{z) + Z 2 (z))(DZ 1 (z) + DZ 2 (z) + zZ'^z) + zZ 2 (z)) 

z 

We solve this linear equation for WZi and WZ 2 and obtain an expression for 
WZ(z) = WZi(z) + WZ 2 (z) in terms of Z x and Z 2 . Then, we replace Z x {z) by 
and Z[(z) by 3Z ^ Z ^ Z ) _ Finally implicit differentiation of the 

functional equation z 2 — zZ 2 (z) + zZ 2 (z) 2 + Z 2 (z) 4 = yields 

2z-Z 2 (z) + Z 2 (z) 2 
21 J z-2zZ 2 (z)-4Z 2 (z)3' 

so WZ can be written in terms of Z 2 and z only. We obtain an expression of 
the form 

N 

WZ(z) 



z3(z 2 - zZ 2 (z) 2 - 3Z 2 (z) 4 ) 2 (z - 2zZ 2 (z) - 4Z 2 (z)3) 2 ' 

where N is a polynomial in Z 2 and z. The denominator only vanishes at z = 
(which is clearly no singularity) and at z = ±13 ^ 3 ~~ 35 , which are singularities 
of Z 2 . Therefore, we only have to consider the expansion of WZ around the 
dominating singularity z = 13v ^|~ 35 , which is given by 

wz ( z)^ 91 

1248 (l - 
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This yields 



(91 - 5\/T3)V^ 5/ 2 ^35 + 13\/l3 X 



E(W n Z n ) ~ ^ » 



312 \ 54 



3.3 p- and Wiener index 

Again, all steps are almost analogous to the previous section. We start from 
the following functional equations, which are simplified by the relations T{z) = 

1-T(z)> K H Z i - l~fi 1 (z)-T(z) and H 2{Z) - (i-T(z)y> • 

zCDRUz) + D(z)) 

= R ^ R f +D ^ + zR [ (z)-R 1 (z), 

DR. W = + + (,) - ft (,) 

_ 2^)T(z)ft(,) + T(,) 2 (DR l( ,) + DR 2 (,)) + ^ _ ^ 



WRi(z) =DRi(z) + 
+ 



z z 
z(WRi(z) + W{z)) 



(l-i? 1 (z)+T(z))2 
2z(zR' 1 (z) + zT'(z))(DR!(z) + D{z) + zR[(z) + zT{z)) 



DRi(*) + 



(i-RM-nzw 

i2i(s) 2 (WRi(*) + W(«)) 



z 

2R 1 {zf{R' l (z) + T , (z))(DR 1 (z) + D{z) + zR[(z) + zT'(z)) 
WR 2 (z) = DR 2 (z) + ^(zmizH^z)) + z(WR l( z) + WR 2 (z)) 



(1-T(z))3 (l-T(z)r 
2z(DRi(z) + DR 2 (z) + zR[(z) + zR 2 (z))zT'(z) 

(l-T(z)f 
2z(zR' 1 (z) + zR' 2 (z))(D(z) + zT'(z)) 

(l-T(z)r 
%z 2 T'{z){D(z) + zT'jz^jR^z) + R 2 {z)) 
(l-T(z)r 



= DRa(z) + 2T(zfW(z)(R 1 (z) + R 2 (z)) + T(znWR 1 (z)+WR 2 (z)) 

z 2 z 
| 2r(z) 3 (DR!(z) + DR 2 (z) + zR'^z) + zR' 2 {z))T (z) 



+ 



2T{zf{R' l (z) + R' 2 {z)){D{z) + zT'(z)) 



z 

6T(z) 4 T'(z)(D(z) + zT'iz^R^z) + R 2 (z)) 
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In this case, we solve the linear equation for WRi and WR2 and insert the 
exact expressions for T, D, W, R\ and R 2 , which can be determined by simple 
quadratic equations: we have 

Ri(z) = J ^1 + Vl - 4i - ^2(1 - 10z + \/T^4i)^ 

and 



R 2 (z)= 1 Vl 4 z -fli(z). 

Note that we always have to take the branch whose value is at z = 0. 
We will use q\ as an abbreviation for \J\ — Az and q 2 as an abbreviation for 



2(1 - 10z + VI -4z). Then we also have 



We obtain an exact expression for WR(z) = WRi(z) + WR 2 (z) in terms of a 
rational function in q x and q 2 . The denominator of this expression is given by 

^(5 gi - 3)(3- 5 gi +g 2 ) 2 , 

which only vanishes at z = 0, z = \ and z = Furthermore, gi has its only 
singularity at j, and q 2 has singularities exactly at j and ^g. Therefore, the 
singularity of WR we have to investigate is 4r . We obtain the expansion 



WR(z) 

which yields 

E(W n a n ) r 



25 ' 
1 



15(1-^) 



VtF 6/2 / 25 



15 V 16 



r(a n ,Z n ) ~ 


(-1.01706) 


• (0.99405)" 


r{cr n ,p n ) ^ 


(1.05088)- 


(0.99023)" 


r(Z n ,p n ) ~ 


(-1.08924) 


• (0.97853)" 


r(W„,cr„) ~ 


(-0.27891) 


• (0.99767)" 




- (0.40351) 


(0.99637)" 


r(W n ,p n ) ~ 


(-1.78357) 


• (0.98209)" 



Table 1: Asymptotic formulas for the correlation coefficients. 
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4 Conclusion 



From the expected values and variances which were calculated in the preceding 
sections, it is possible now to deduce the asymptotic correlation coefficients 
for the investigated indices (Table ^) . For their interpretations and further 
discussion see [TT] . 
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